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Fermi-LAT reveals two huge gamma-ray bubbles existing in the Galactic Center, called ’Fermi 
Bubbles’. The existence of two microwave bubbles at the same region are also reported by the 
observation by WMAP, dubbed ’WMAP haze’. In order to explain these components. It has been 
argued that the gamma-rays arise from Inverse-Compton scattering of relativistic electrons accel¬ 
erated by plasma turbulence, and the microwaves are radiated by synchrotron radiation. But no 
previous research reproduces both the Fermi Bubbles and WMAP haze under typical magnetic fields 
in the galaxy. We assume that shocks present in the bubbles and the efficiency of the acceleration 
by plasma turbulence, ’stochastic acceleration’, changes with the distance from the shock front. 

The distance from the shock front increases with time, accordingly the efficiency of the acceleration 
changes with time. We also consider the time development of the electrons escape from the tur¬ 
bulence by diffusive loss. Our model succeed to reproduce both the observed characteristics of the 
Fermi Bubbles and WMAP haze under typical magnetic fields. 


I. INTRODUCTION 


Gamma ray data from the Fermi-LAT reveal 
two giant bubbles extending up to ~ 50° above 
and below the Galactic disk, called ’the Fermi 
Bubbles(FBs)’[l|[3. The spectra of the FBs are 
harder than ambient radiation fields, and they have 
’sharp edge’, which means sudden change of the 
brightness at the boundary. The surface brightness 
is almost constant in the FBs, though it doesn’t mean 
the FBs have constant volume emissivities. There 
also exists the microwave bubbles in the same region, 
dubbed ’the WMAP haze’Q. These huge structures 
may suggest past large-scaled activities in the Galac¬ 
tic Center(GC), near the massive black hole Sgr A*. 

Two representative models have been argued in or¬ 
der to explain both the FBs and the WMAP haze. 
’Hadronic model’ reproduces the FBs by the gamma 
ray from the 7r°-decay process, and the WMAP haze 
by the radiation from the secondary electrons QQ. In 
the Hadronic model, protons are accelerated by the 
Diffusive Shock Acceleration(DSA) or other process, 
and they make after collision with cosmic-ray pro¬ 
tons. TT^ are also made in this process, and they decay 
to e^. These secondary particles radiate microwave in 
the magnetic field by the synchrotron radiation. 

In contrast, ’Leptonic model’ explains the FBs by 
the gamma ray from the Inverse Compton scatter- 
ing(IC) process, and the WMAP haze by the sto- 
chrotron radiation from the relativistic electrons QQ- 
In the Leptonic model, it is assumed that electrons 
are accelerated by the second order Fermi accelera- 
tion(stochastic acceleration) by plasma turbulence or 
molecular clouds, in order to explain hard spectra of 
the FBs. These relativistic electrons give their en¬ 
ergy to the ambient photons by the IC process, and 
they are observed as gamma ray. Some electrons radi¬ 
ate microwaves by the synchrotron radiation, and the 
WMAP haze are made. 

However, both Hadronic and Leptonic model have 


problems. Hadronic models require so much energy ~ 
10®^ erg, and estimated microwave spectrum conflicts 
with the WMAP dataQ. In the leptonic models, it 
is difficult to explain the WMAP haze under typical 
galactic magnetic field(i? ^ 4^G), so it is required 
strong magnetic field in the large scale Q- 

So we construct the extended model based on the 
stochastic acceleration model by Mertsch & SarkarQ. 
Their model assumes that electrons are accelerated by 
the turbulence whose acceleration efficiency decreases 
with the distance from the shock front. But they don’t 
consider the time development of the electron spec¬ 
trum, and presume that electrons are at steady states. 
They also neglect the contribution of the electrons es¬ 
cape from the turbulence. So we extend their model 
by considering both time development and electrons 
out of the turbulence. In the secHIl we explain our 
acceleration model. And in the sec lllli we show the 
results of our numerical calculation. 


II. MODELS 


We assume thet there is the shock front near the 
edge of the FBs and the plasma turbulences are drifted 
from the shock front with time by the advection at 
V = Vpro- 

In order to calculate the electron spectrum, we solve 
the Fokker-Planck equation, 


dn d f 2 j-, d n\ n d 
dt dp V dpp^)^ tesc dp 

where the momentum diffusion 
the second term is described Q 



Qinj — 0 

( 1 ) 
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which represents the farther from the shock, the 
weaker the acceleration efficiency becomes. 

In ([2]) p is energy, k is wavenumber, vf{^) is 
the velocity of the fast-mode wave, and W{k,^) = 
(MV47r)L-2/3fc-ii/3 is the energy density of the Kol¬ 
mogorov turbulence. DxxH) is the spatial diffusion 
coefficient. The dimensionless parameter ^ = xjLix 
is the distance from the shock front, and L is the tur¬ 
bulence size) represents how far from the shock front. 

The third term describes the diffusive loss from 
the turbulence, and then t^sc is accounted for t^^^c = 

LyOxx- 

The forth term of ([1]) corresponds to the energy loss 
of the synchrotron radiation and the IC scatering. 


( ^ 
\dt. 

\ — Psyn + PlC 

‘ cool 

(3) 
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Pic = 

-ca^yyUph X /kn 

(5) 


Ub = B^jSn where i? = 4 jiG, and Upt is shown 
in Figured] /kn < 1 is the factor due to the Klein- 
Nishina effect,which is numerically evaluated. In our 
models, fpro is fixed Vpro = 250kmsec“^, which cor¬ 
responds to sound speed at kT ^ keV. 

We presume that injection rate is constant, 

Qinjip) = QoS{p-po) ( 6 ) 


where the spectrum shape of injected electron is as¬ 
sumed to be ^-function, ninj{p) (x 5{p — po)- In this 
paper,we use cpo = 10® eV. 



FIG. 1: Background Photon Fields^.We use the black 
solid lme(i? = 2 kpc ,2 = 5 kpc) as Uph- 


We consider the effect of time development, so we 
assume that the distance from the shock front changes 


with time like ^ = xjL = Vprot/L, and then Dpp be¬ 
comes the function of the time t. We also consider the 
contribution of the electrons out of the turbulence, so 
we solve the equation of these electrons 


duesc 9 (dp _ ^ 

^ + Wp 




(7) 


where is equal to the escape rate from the tur¬ 
bulence = n/tesc- 

To verify the effect of the escape, we calculate the 
spectrum of electrons under three cases. ’No escape 
model’ is the model in which we neglect the escape 
term n/tesc(the limit t^-sc —t oo). In the ’Cut escape 
model’, we consider the effect of escape but we don’t 
deal with electron out of the turbulence(like[3|). So 
in the Cut escape model, we solve only o, and don’t 
solve 0. We consider the electrons out of the tur¬ 
bulence in the ’With escape model’, then we solve 
both o and 0, and finally calculate sum of them 

rie — n Tlesc- 

After the calculation of the spectrum of the elec¬ 
trons, we estimate the Intensity from those electrons 
and compare with the data of the Fermi-LAT and the 
WMAP. In the next section sec lIIIl we show the re¬ 
sults of our calculation. 


III. RESULTS 

Our results are shown in the Figure [D We cal¬ 
culate equation o under the three cases at variable 
parameters L and Qq for t = 2kpc/r'pro — 8 Myr. 
Vpro = 250kmsec“^ and B = 4/rG are fixed. The result 
of ’No escape model’ is shown in blue solid line, ’Cut 
escape model’ is shown in red solid line, and ’With 
escape model’ is shown in black solid line. The data 
points in are adopted. 

Gamma ray data of the Fermi-LAT are nicely re¬ 
produced in all cases. But the microwave data of 
the WMAP are reproduced only by ’With escape 
model’(Figure m. This is because of increasing of low 
energy particles due to escape from the turbulence. In 
this case the turbulence size L = 10 pc, smaller than 
that of in 0 L = 2 kpc. If we reproduce the WMAP 
data in other cases, stronger magnetic field is required 
at large scale. 


IV. SUMMARY 

We reproduce both the Fermi Bubble and the 
WMAP haze under typical galactic magnetic field 
B = ApG by considering the effect of time dependence 
and electrons which escape from the turbulence. ’Cut 
escape model’ which is similar to the model of previ¬ 
ous research Mertsch & Sarkar Q reproduce gamma 
ray data well, but predicted microwave emission is not 
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Photon Spectrum Photon Spectrum 



FIG. 2: left: Photon spectrum under three cases, right: Photon spectrum at gamma ray region. ’No escape moder(blue 
solid), ’Cut escape moder(red solid), ’With escape moder(black solid) are all shown.The data points are from [^. 


enough at B = 4/rG. Our results suggest that smaller Acknowledgments 
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